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Abstract. We theoretically demonstrate a counter-intuitive phenomenon in optical 
interferometry with a thermal source: the emergence of second-order interference 
between two pairs of correlated optical paths even if the time delay imprinted by 
each path in one pair with respect to each path in the other pair is much larger than 
the source coherence time. This fundamental effect could be useful for experimental 
simulations of small-scale quantum circuits and of 100%-visibility correlations typical of 
entangled states of a large number of qubits, with possible applications in high-precision 
metrology and imaging. As an example, we demonstrate the polarization-encoded 
simulation of the operation of the quantum logic gate known as controlled-NOT gate. 
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1 . Motivation 

The Hanbury Brown and Twiss (HBT) effect [1, 2] discovered in 1956 was at the heart of 
the development of the held of quantum optics. Indeed, this discovery led to numerous 
remarkable multiphoton experiments which have not only deepened our fundamental 
understanding of multiphoton interference [3, 4, 5, 6, 7, 8, 9] but have also led to 
numerous applications in information processing [10, 11, 12, 13, 14, 15] and imaging 
[16, 17, 18, 19, 20, 21, 22], 

The HBT effect fundamentally reveals the second-order coherence of a thermal 
source. For example, second-order temporal correlations can be measured after the 
interaction of multi-mode thermal radiation of given bandwidth Au with a beam splitter: 
two detectors at the beam splitter output ports have twice the chance to be triggered 
at equal times than with a relative time delay longer than the coherence time XjAbj of 
the thermal held. 

It is interesting to modify the described HBT scheme by adding two unbalanced 
Mach-Zehnder interferometers at the beam splitter output ports, the “control” port 
C and the “target” T, as depicted in hgure 1(a). We further consider path lengths 
L ■■= Lc = Lt and S ■= Sq = St such that the time | L — S' | /c, with c the speed 
of light, is much larger than the coherence time of the source. Can we observe second 
order interference by performing correlation measurements at equal detection times at 
the output of the two Mach-Zehnder interferometers? 

Interestingly, in this paper we show that second-order interference between the two 
pairs (Lc, Lt) and {Sc-, St) of optical paths {multipath correlation interference) can be 
observed even if the time delay \ L — S \ /c imprinted by each path in the pair (Lc, Lt) 
with respect to each path in the pair {Sc, St) is much beyond the source coherence 
time. Furthermore, we describe the similarities and differences between this second- 
order interference effect and the one demonstrated by Franson [7] in 1989 by using a 
two-photon entangled source. 

The fundamental interference phenomenon described here is also of interest in view 
of the recent studies of simulations of quantum logic operations and entanglement 
correlations using classical light [23, 24, 25, 26, 27, 28, 29, 30]. In particular, 
by considering the interferometric scheme in figure 1(b), we demonstrate how this 
interference effect is able to simulate the result of a controlled-NOT (CNOT) logic 
operation [10, 31, 32, 33, 34, 35, 36, 37]. 

The paper is outlined as following: we demonstrate how multipath correlation 
interference can be observed with a thermal source in section 2; we apply this novel 
phenomenon to the simulation of a CNOT gate operation in section 3; we show how this 
second-order interference effect can be generalized to interferometers based on arbitrary- 
order correlation measurements in section 4; and we conclude with discussions in section 
5. 
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Figure 1. (a) Modified HBT interferometer with a thermal source where two 

unbalanced Mach-Zehnder interferometers are placed at the output channel of the first 
beam splitter and the correlation between the photon-number fluctuations is measured 
at the interferometer output, (b) Interferometer simulating a CNOT-gate operation: 
1) The first transformation Wprep prepares the initial polarization state; 2) By using two 
unbalanced interferometers a polarization-dependent transformation P is implemented; 
3) At the interferometer output the polarization-dependent correlation between the 
photon-number fluctuations is measured. 


2. Multipath Correlation Interference 


We consider here the interferometer in hgure 1(a). The interferometer has only one 
sonrce, which generates in the input port A thermal light with a given horizontal 
polarization H. Therefore, the input state is described by [38, 39] 
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with the Glauber-Sudarshan probability distribution [3, 40] 
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where is the average photon number at frequency oo. For simplicity, but without 
losing generality, we consider a Gaussian frequency distribution [38] 
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with the mean photon rate r, average frequency cuq and spectral width Aoj. 

At the interferometer output correlation measurements in the photon-number 
fluctuations 
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at the detection time t^ around the mean value (nf), with d = C,T, are performed by 
using either photon-number resolving detectors or single-photon detectors [41, 42] with 
integration time 6t l/Zlo;, which are currently available (e.g. ~ 10 — 100ns). The 

expectation value for the product of the photon-number fluctuations in (2) at the two 
output ports is [41, 42, 38] 

{AucAut) = {nc ut) - {ric) {ut) 

(X G^‘^\tcOT) — G^^\'tc^tc)G^^\'tTOT) = \G^^\'tcOT)\ ■ (3) 

Here we used the properties [38] 

{ucUt) oc G^‘^\tc,tT), 


with 

G^‘^\tc,tT) = G^^\tc,tc)G^^\tT,tT) + \G^^\tc,tT)\ (4) 


and 

(Ud) oc G^^^tdCd), 


where and are, respectively, the first and second-order correlation functions 
[38]. Therefore, the outcome (3) of the correlation measurement does not depend 
on the “background” term G^^\tc,tc)G^^'^{tT,tT) in (4). As we will show, the only 
relevant term \G^^\tc,tT)\ characterizes the second-order interference occurring in the 
interferometer. Toward this end, we introduce the definition of the hrst-order correlation 
function 


tr PhE^c \^c)E^\tT) 


(5) 


Here, in the narrow bandwidth approximation. 



( 6 ) 


is the electric field operator at the detector d = G, T in terms of the frequency-dependent 
annihilation operator dA{oj) at the only input port A where a source is placed; while 
E^\td) is its respective Hermitian conjugate. The factors ^nd describe 

the propagation through the paths Ld and Sd, respectively. By using (5) and (6), and 
defining the effective detection times 
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with Id = Ld, Sd, we show in Appendix A that the expectation value in (3) can be 
written as 
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with a constant, describing the contribntion of the corresponding path pair {Ici^t), 
depends on the time delay between the detected photons in the 

two output ports before the propagation through the two paths Ic and It, 
respectively. In a standard HBT experiment only a fixed pair of paths can 
contribute to the measurement. Differently, here all the four possible pairs {Ic, It) = 
{Lc, Lt), {Sc, St), {Lc, St), {Sc, Lt) of paths can lead to a joint detection and, in 
principle, can interfere as in (8). 

Now, we consider the following conditions for the delays between the effective 
detection times defined in (7): 
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which, by using (10), corresponds to differences between the path lengths in each Mach- 
Zehnder interferometer much larger than the coherence length cjAuj of the source: 


Lc — iS'c'l 3> c/Au, \Lt — 5't| ^ c/Au. 


( 12 ) 


The conditions (10) can be simply achieved experimentally, for example, in the limit 

tc — tn ^ \tc ~ ^r| ^ 1/Au (13) 

of approximately equal detection times with respect to the coherence time XjAbj and 
approximately equal paths 


T ’■— Tc — Tt aa \Lc — Lt\ ^ c/Au!, S ■ — Sc — St I'S'c — 'S'j'l c/Au 

(14) 

with respect to to the coherence length c/Au. By using (14), the two conditions (12) 
reduce to the single condition 


|L — S'! S> c/Au. (15) 

In the conditions (11) and (10) and by using (9) the expression in (8) becomes 

{AucAht) oc |G(^^(tc,^r)f = ^t) + 

= ‘la^r^ (1 + cos ipL-s), (16) 

with the relative phase 

[(tA'’ - M’) - (tC - *7"’)] = T [(io - Lt) - {Sc - St)1 . (17) 


where we used (7) in the second equality. The expectation value (16) depends now 
on the interference between only two contributions and associated 












Multipath Correlation Interference and Controlled-NOT Cate Simulationwith a Thermal SourceG 

with the pairs of optical paths {Lc,Lj’) and {Sc, St), respectively. The paths Lq and 
Sc in the interferometer are correlated with the paths Lt and St, respectively, and 
only the corresponding path pairs {Lc, Lt) and {Sc, St) interfere. This interference 
occurs even if the differences \Ld — Sd\, with d = C,T, between the path lengths in 
each Mach-Zehnder interferometer are much larger than the coherence length of the 
source (see (12)). Indeed, for a thermal source, the interfering contributions 
and (see (9)) do not depend on the relative path lengths in each Mach-Zehnder 

interferometer but only on the difference between the delays and 

between the detected photons in the two output ports before the propagation through 
the two pairs {Lc, Lt) and {Sc, St) of paths, respectively. Since these time delays are 
very small compared to the coherence time of the source (see (10)) both pairs {Lc, Lt) 
and {Sc, St) contribute to the observed second-order interference. 

It is worthwhile to compare the interferometer described here with the famous 
Franson interferometer [7], where the state at the output of the beam splitter in 
figure 1(a) is substituted by a two-photon entangled state and the coincidence rate 
for detecting a photon at the same time in both output ports is measured. For a full 
comparison, we hrst determine the coincidence rate associated with the absorption of 
a single photon from the field at each of the two output ports of the interferometer 
in hgure 1(a) considered here. This corresponds to measure the standard second-order 
correlation function [38] in (4) at approximately equal detection times tc = tT (see (13)) 
and in the conditions (14) and (15) for the interferometric optical paths. In particular, by 
adding the product G^^'^{tc,tc)G^^'>{tT,tT) = 4a^r^ of the intensities at the two output 
ports to the second-order interference term |G*^^^(fc',t'r)| found in (16), we obtain a 
second-order correlation function 

G^‘^\tc, tT) OC 3 -|- cos (fL-S 

with visibility 1/3. This is a crucial difference between the interferometer described here 
and the Franson interferometer, where the second-order correlation function measured 
at the output at equal detection times manifests second-order interference with 100% 
visibility. Instead, in the interferometer in figure 1(a), 100%-visibility interference is 
only achieved by measuring the correlation (16) between the fluctuations in the number 
of detected photons, where the “background” constant term G^^\tc,tc)G^^\tT,tT) is 
effectively “subtracted” from the second order correlation function. Therefore, the 
emergence of this interference effect is very different from the physics of two-photon 
interference based on energy-time entanglement in the Franson interferometer. Indeed, 
in the Franson interferometer the interference between the two pairs {L, L) and {S, S) 
of optical paths emerge from the fact that the two input entangled photons are emitted 
at the same time and the joint emission time is uncertain in the quantum sense. 
Therefore, 100%-visibility second-order interference can be observed even if the first- 
order coherence length is much less than the difference L — S between the path lengths. 

In the interferometer described here, instead, neither an entangled source is used nor 
any entanglement process occurs. Therefore, coincidence events do not necessarily 
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correspond to the absorption from the field of photons which entered the interferometer 
at the same time. In principle, the detected photons could have taken any of the 
four possible pairs of paths from the source to the two detectors. Nonetheless, the 
interference term \G^^\tc,tT)\ in the second-order correlation function, emerging from 
the measurement of the correlation (16) between the fluctuations in the number of 
detected photons as an average over all the possible experimental outcomes, contains 
only the contribution of two indistinguishable pairs {Lc, Lt) and {Sc-, St) of correlated 
paths |. In particular, the interference pattern (16) depends on the difference 

Tl-s = Tc — Tt 

in (17) between the relative phases ipd = oJo{Ld — Sd)/c, with d = C,T, in the two 
Mach-Zehnder interferometers. Differently, in the Franson interferometer the resulting 
interference pattern depends on the sum of the relative phases in the two Mach-Zehnder 
interferometers. 


3. CNOT Gate Simulation 


The interference phenomenon based on multipath correlations demonstrated in the 
previous section can be used to reproduce on-demand correlations in different degrees of 
freedom without the use of entanglement. Here, we address, for example, the simulation 
of a CNOT gate operation by encoding these multipath correlations in the polarization 
degree of freedom. 

For this purpose, we consider the interferometer in figure 1(b). The source at the 
input port A is again described by (1). The i7-polarized thermal light impinges on the 
balanced beam splitter BSl and, by using half-wave plates, is prepared at the “control” 
port C and “target” port T in two general polarizations (pQ = {cos pc smpc) and 
<Pt = (cos 07- sin 07-)^, respectively. Here, the H and V polarization directions are 
indicated by the vectors (1 0)^ and (0 1)^, respectively. Thereby, by introducing the 

polarization rotations 

cos 0c,r sin 0c,r 

^sin 0c,r - cos 0c,r^ 

the interferometer transformation associated with the first part of the interferometer 
connecting the input ports A, B with the ports C, T is given by 


R 


4>c,t 


^prep ■ 


^ 1 7 y — ^ I 



(18) 


where we used the expression for the balanced beam-splitter transformation Ubs, with 
1 := diag(l,l). 


f Interestingly, there is, in principle, no upper bound to the difference \L — S\ between the path lengths 
in (15) which limit this interference effect for an ideal stationary thermal source. This is of course not 
the case in “real world” experiments. 
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The second part of the interferometer consists of a “control” interferometer 
connecting the ports C and C and a “target” interferometer connecting the ports T and 
T. Therefore, the global interferometric evolution is described by the two polarization- 
dependent transformations Vc,c ^md Vt,t dehning the diagonal matrix 

V ■■= diag {Vc,c, ■ (19) 


In particular, in the control interferometer, the light in the polarization modes H and V 
at the output of the hrst polarizing beam splitter PBSl acquires the time delays Sc/c 
and Lc/c, respectively, with c the speed of light, before being recombined at the output 
of the second polarizing beam splitter PBS2. This leads to the control transformation 

Vc,c ■■= diag {C‘^Sc/c^^iu.La/cy 


On the other hand, the light in the target interferometer is coherently split into 
two different paths and recombined by the balanced beam splitters BS2 and BS3, 
respectively, independently of the polarization. The light polarization is unchanged in 
the path of length St- Instead, in the path of length Lt the polarization modes H and 
V are flipped (IT -H- V) by the NOT-gate operation F ■■= antidiag (1,1) implemented 
by a half-wave plate with axes rotated by vr/d with respect to the H and V axes. Thus 
the overall target evolution is described by the transformation 

Vt,t ■■= ^ . 

By using (18) and (19), we derive the total interferometer matrix 


M(uj) 


I _ f i T^C,C R<l)e Rc,C R<j)c \ 

\/2 yPr,r-R(^r ^Rt,t R<f>r) 


( 20 ) 


We dually address the detection process, consisting of measuring the polarization- 
dependent correlation between the fluctuations in the number of photons Ang^(tc) 
and Anej,(tT) detected at the control and target ports d = C,T, respectively, with 
polarization 6^ ■= (cos 9d sin 9d) at time td- The expectation value for the product 
of the photon-number fluctuations at the two output ports is [41, 42, 38] 


(Ane^Ane^) oc 
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where the hrst-order correlation function 
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is now, differently from (5), polarization dependent. Here, in the narrow bandwidth 
approximation and in the (H, H)-polarization basis. 
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with the elements ^ = M.c,Ai M.t,a in the first column of the total interferometer 
matrix (20), is the electric field operator at the detector d = C, T in terms of the 
frequency-dependent annihilation operators (a;) and (cu) at the only input port 
A where the source is placed, while (td) is its respective Hermitian conjugate. In 
Appendix B we show, that, in the limits (11) and (10), (21) becomes 


{Ang^Ane^) oc {tcOT) 

oc I cos (pc cos 9c cos ( 07 — 9t) + sin (pc sin 9 c sin {<pj- 9 t) I 


{Sc, St)/ 


( 22 ) 


with ifL-s defined in (17). 

Here, similarly to a CNOT gate operation, the control path Sc associated with the 
polarization mode H is correlated with the target path St where the light polarization 
remains unchanged; instead, the control path Lc associated with the polarization mode 
V is correlated with the target path Lt where a NOT-gate operation occurs. Moreover, 
analogously to the scheme in figure 1 (a), the two pairs (Lc, Lt) and (Lc, St) of optical 
paths interfere. 

Let us fully compare now the interferometer in hgure 1(b) with a genuine CNOT 
entangling operation on the two-qubit input state | 0 c)c \4>t)t^ where 


\(pc)c ■■= cos 0 c \H)c + sin 0 c \V)c , 


and 

I0r)r cos0r |Lf)^ + sin0r |H).^ 

are expressed as superpositions of the polarization states |Lr) and \ V), corresponding to 
single-photon occupations of the H and V modes, respectively. A CNOT-gate operation 
on this input state leads to the output entangled state 

{F)\ 


)cT = COS (pc |iL)c |0r)T + sin0c \V)c Wr )t^ 


where 


( P'] 

\(p}j- 0^ ■.= sin 07 - \H)j, -|- cos 0 r \ V)rp . 

'>C,T 


Polarization-correlation measurements over the state 10)^^- occur with a probability 


P( 


CNOT •- 


(6'c)^t|0)ct ~ |cos0c cos 6*c cos ( 07 — 6 't) + sin0csin6*csin (07--|-6 't)|^ . 

(23) 


Comparing (23) with (22) in the limit (pi-s 1, we obtain the expectation value 


(Ang^Ang^) oc r^PcNOx, 


(24) 


which takes into account all the possible outcomes for the product of the photon-number 
fluctuations measured at the output of the interferometer in hgure 1(b). We emphasize 
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Figure 2. iV-order interference network with a thermal source based on: 1) 
Preparation of the A^-channel initial state for arbitrary polarization angles (pi [i = 
1,... ,iV) with a generalized transformation Uprep implemented by a symmetric 2N- 
port beam-splitter and N half-wave plates; 2) A^-channel polarization-dependent 
evolution consisting of polarization rotations and P-type transformations as in 
figure 1 (b); 3) Polarization-dependent correlation measurements in the photon-number 
fluctuations at the output of the network. 


that no entanglement process occnrs in the interferometer; therefore the proposed 
scheme is not an entangling gate. Nonetheless, the measnrenient of the correlation (24) 
between the photon-nnmber fluctuations at the two output ports allow us to simulate a 
CNOT-gate operation. As a “bonus”, we find that the correlation signal can be enhanced 
on demand by increasing the square of the source mean-photon rate, making it robust 
against technological losses. 

As an example, if we fix the polarization angles (fc = 7r/4 and fj- = 0 in the setup 
in hgure 1(b), the expectation value in (24) reads 

(Ane^Ane^) oc cos^ (9c - 6 t) , 

simulating the 100%-visibility correlations typical of a Bell state = 

(\HH) + \VV))/V2 even if no Bell state is produced. Indeed, neither an entangled 
source is used, as for example in the experiment of Sanaka, Kawahara and Kuga 
[37], nor an entanglement process occurs in the interferometer as in a genuine CNOT 
gate operation. Nonetheless, an observer T can provide to a separate observer C the 
measured fluctuations Ang^ in the number of photons for a polarization angle 6c (in 
a given computational base) unknown to T, who can infer the value of 6c only based 
on his/her corresponding measurements of Ang^ (with 6t in the same computational 
base). Interestingly, differently from a genuine Bell state, the correlation between the 
measurements of the two observers C and T emerge only from the expectation value of 
the product of the corresponding photon number fluctuations Ang^ and Ang^.. 

4. N-order Interference Networks 

One can finally generalize the second-order interference effect described in this paper 
to arbitrary orders N. For example, we consider an A^-order interferometer as in 
figure 2 which generalizes the scheme in hgure 1(b). For this purpose, we hrst prepare 
the A^-channel initial polarization state: i7-polarized thermal light impinges on a 
symmetric 2A^-port beam-splitter (generalization of the balanced beam splitter BSl 
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in figure 1(b)); at each output port the polarization is then rotated, respectively, 
by arbitrary polarization angles 0*, with i = 1,2, by using N half-wave 

plates. This prepared initial state propagates through a 2A^-port polarization-dependent 
interferometer, consisting of polarization rotations and P-type transformations as in 
figure 1 (b). Finally, the polarization-dependent correlation between the photon-number 
fluctuations is measured at the N output ports at approximately equal detection times. 

We emphasize that for N = 2 one may also measure this correlation (24) by 
subtracting the product of the independent light intensities at the two detectors from the 
standard second-order correlation function. Differently, in higher order interferometric 
networks it may be useful to directly measure the correlation between the photon- 
number fluctuations. 

Experimentally, the precision in measuring the expectation value of the product 
of the photon-number fluctuations Ang^ at each output port d = 1,2, ...,N of a given 
network depends on the number n of performed measurements. For large numbers 
n of measurements the distribution of the measured mean values is normal around 
the expectation value with an indetermination given by the indetermination in the 
distribution of the measurement outcomes normalized by the root of n. Evidently, the 
indetermination in the distribution of the measurement outcomes depends intrinsically 
on the given iV-order interference network. In general, a scaling in the number of 
resources typical of a quantum network with genuine entanglement cannot be achieved 
by the scheme described here. Nonetheless, an iV-ordered interferometer as in figure 2 
could be used to simulate experimentally 100%-visibility correlations typical of entangled 
states of N qubits [30], such as GHZ states, as well as small-scale quantum circuits and 
algorithms. 

5. Discussion 

We demonstrated a novel interference phenomenon emerging from the fundamental 
nature of multipath correlations with a thermal source. In particular, we introduced a 
novel interferometer (see figure 1(a)) where full correlations between the interferometric 
paths Sc and St {Lc and Lt ) emerge at the interferometer output from measuring the 
correlation between the fluctuations in the number of detected photons at the two output 
ports. We showed how the interference between the two pairs {Sc, St) and {Lc, Lt) of 
correlated optical paths occurs even if the time delay imprinted by each path in one pair 
with respect to each path in the other pair is much larger than the source coherence time. 

We also pointed out the differences and the similarities between this interferometer and 
the well known Franson interferometer where an entangled two-photon source is used 
instead of a thermal source. 

The interference effect demonstrated here can be easily observed experimentally: 1) 

It relies on one of the most natural sources which can be easily simulated in a laboratory 
by using laser light impinging on a fast rotating ground glass [43]; 2) The calibration 
of the interferometric paths and of the detection times (see (15), (14) and (13)) can 
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be easily achieved in the case of a thermal source, where the coherence time can range 
from the order of ns to fis. In an analogous way, multipath correlation interference can 
be also experimentally observed in the spatial domain with a spatial-mode dependent 
thermal source [44], 

In conclusion, this interesting phenomenon provides a deeper fundamental 
understanding of the physics of coherence, multipath correlations and interference using 
a thermal source. Furthermore, it can be used to implement correlations in different 
degrees of freedom without recurring to entanglement processes. As an example, we 
demonstrated the polarization-encoded simulation of a CNOT logic operation. We also 
showed how this second-order interference effect can be extended to arbitrary orders 
N, leading to the interference of more general conhgurations of correlated paths. This 
could be used to simulate on-demand 100%-visibility correlations typical of entangled 
states of N qubits, with possible applications in high precision metrology and imaging 
[19, 21, 45, 22], and in the development of novel optical algorithms [46, 47, 48, 49, 50]. 

Acknowledgments 

V.T. would like to thank M. Cassano, M. D’Angelo, P. Facchi, J. Fan, J. Franson, M. 
Freyberger, A. Garuccio, Y.-H. Kim, S. Laibacher, A. Migdall, S. Pascazio, T. Peng, 

T. Pittman, W.P. Schleich, and Y.H. Shih for fruitful conversations about the proposed 
scheme. V.T. is also grateful to M. Cassano, T. Peng and Y.H. Shih for discussions, 
during his visit at UMBC in 2014, about the experimental realization of this proposal 
in the spatial domain in Y.H. Shih’s laboratory (paper in preparation). 

V.T. acknowledges the support of the German Space Agency DLR with funds 
provided by the Federal Ministry of Economics and Technology (BMWi) under grant 
no. DLR 50 WM 1136. 


Appendix A. Correlation between the photon-number fluctuations for the 
interferometer in figure 1(a) 


We And here the explicit expression of the expectation value (in (3)) 

{AncAnr) (x\G^^\tc,tT)f (A.l) 


of the product of the photon-number fluctuations by calculating the hrst-order 
correlation function (in (5)) 




tr 


PhE^ \tc)Ej^XtT) 


(A.2) 


with the electric held operators 


E^itc) 



(A.3) 




Multipath Correlation Interference and Controlled-NOT Cate Simulationwith a Thermal Sourcel?) 


and 




(A.4) 


where is a constant. By inserting (A.3) and (A.4) into (A.2), and using the dehnition 
(in (7)) of the effective detection times 


Aid) _ ^ 


we obtain 




iA:^ 


tr ( Pi, /dw’*( 1 ^) 0 ^"^' 


+ tr ( PH [dLj fdui’e'^A’ A) 


+ tr I pHfdu,fdu,'d(A^°'--A-\A’\u,)af\j) 


+ tr I pHfdu,fdu,'d(A^°'--A-\Apu,)a'-"\u,') 


By using the property [9, 38] 

for the multimode thermal state pn in (1), we obtain 




G^^\tcOT) = - 




0"‘~UL‘'C J0 \-^C ^ 0 ‘'T J0 

which can be rewritten as 


with the contributions 


G^^\tcOT) — G^^^’’‘'^\tc, tx), 

IcjIt 


(A.5) 


for all possible pairs (IcJt) = {Lc, Lt), {Sc, St), {Lq, St), {Sc, Lt) of optical paths. 
By substituting (A.5) in (A.l) we hnally hnd the expression of (AncA^UT) in (8). 

Appendix B. Polarization-dependent correlation between the 
photon-number fluctuations for the interferometer in figure 1(b) 

We calculate here the expectation value (in (21)) 


{AnecAne^)c^T^ G^Pcfir^'^c^G) 


(B.l) 
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for the product of the photon-number fluctuations, with the polarization-dependent 
hrst-order correlation function 


^(1) 

^8 c ,0 t 


(tcytr) 


tr 


h(-) 


Ph [&c ■ (tc)) [Ot ■ Erp (t 


^(+)/ 


(B.2) 


The electric held operator E^^\tci), with d = C,T, for direction of propagation 
perpendicular to the H-V plane is given in the narrow bandwidth approximation by 
the operator 

with a constant K, the elements Aid,A = Me, a, Mt,a in the hrst column of the total 
interferometer matrix A4(cj) in (20), and the annihilation operators and a^^(ci;) 

at the only port A where a source is placed, while \'ld) is its respective Hermitian 
conjugate. 

By dehning 


Ac ■= (cos 9c sin 9c) M.c,a{^) ( g j ^ + sin 9c sin , 


and 

At ■= {cos 9 t sin 6 *t) M.t,a{^) 
1 


0 


2^2 


[cos^T (cos+ sin^T (sin, 


it is useful to introduce the effective held operators 


£^/\td) ■■= iKldu Ad e-^‘^^‘^d^^>{u), 




(B.4) 


with d = C,T. Indeed, given the hxed polarization H of the thermal light produced by 
the source in the port A, equation (B.2) can be rewritten by using (B.3) and (B.4) as 


^(1) 

^dcdT 


{tcAr) 


tr 


Ph (^c \'tc) 



By direct substitution and using again the dehnition (in (7)) of the ehective detection 
times 


Aid) _ y. 

td -trf ^ , 


(B.5) 
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with Id = Ld, Sd, we obtain 




iK^ 


cos 9c cos fc (cos 9t cos fr + sin 9t sin fj) 

X tr i^HJduj (cn')^ 

- cos 9c cos 0 c (cos 9t sin 07 - + sin 9t cos 0 r) 

X tr 

- sin 9c sin 0 c (cos 9t cos 0 r + sin 9t sin 07 -) 

X tr (p„ydt.ydt.'e><“‘?°’-“'‘?"’>a^«>*(x,)a!f'(at')) 

- sin 9c sin 0 c (cos 9t sin fj- -|- sin 9t cos 0 r) 


By using again the property [9, 38] 

tr fp 77 ) = fC-^o{h-t 2 )^-{h-t 2 )^Au^y 2 


for the multimode thermal state pn in ( 1 ), we obtain 
= - -^r 


cos 9c cos (pc cos {9t — 0 r) ^ 
-I- cos 9c cos Pc sin {9t + pr) 




6'c'sin0ccos(6'T - 


+ sin 9c sm pc cos [9t — Pt) c 

-h sin 9c sin pc sin {9t + Pt) _ 

(B. 6 ) 


Finally, by applying the conditions (in ( 11 ) and (10)) 

tg'c) _ pSr) ^ _ pSc) 


ALc) _ ALt) 

Lc irp 


Au S> 1, 
Au <C 1, 


psc) _ ASt) 
he (t 


Au S> 1, 
Au 1 , 


(B. 6 ) reduces to 




_ /^AcAt)(± ± \ , /^{ScAt) X 'i 

~ ^9c,9t [tcAT), 


with the two contributions 

cos Pc cos 9c cos {pr - 9 t) e'‘^° ], 

(^c-, tr) = sin pc sin 9c sin {pr + 9 t) e'“° ]. 
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Thereby, (B.l) reads 


{AnecI^ne^)c^T ^ Gf^l)i'^\tc, tr) + tr) 

which, by introducing the relative phase (in (17)) 

Tl-s ■= [{Lc — Lt) — {Sc — 5'r)], 


finally reduces to 

{AnecAn 0 ^)(jj, oc |cos (fc cos 9c cos (07— 9t) + sin 0c sin 9c sin (07- + 6 *t) , 

as in (22). 
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